Universality of the 1/3 shot-noise suppression factor in nondegenerate diffusive 

conductors 
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Shot-noise suppression is investigated in nondegenerate diffusive conductors by means of an ensem- 
ble Monte Carlo simulator. The universal 1/3 suppression value is obtained when transport occurs 
under elastic collision regime provided the following conditions are satisfied: (i) The applied voltage 
is much larger than the thermal value; (ii) the length of the device is much longer than both the 
elastic mean free path and the Debye length. By fully suppressing carrier-number fluctuations, long 
range Coulomb interaction is essential to obtain the 1/3 value in the low-frequency limit. 



PACS numbers: 72.70.+m, 72.20.-i, 73.50.Td, 73.23.Ad 

In recent years kinetic phenomena in mesoscopic struc- 
tures are offering a fascinating scenario for fundamental 
research jjj. One of the most up-to-date subjects is shot- 
noise suppression in disordered conductors. Here, the ex- 
cess noise power has been predicted to comprise exactly 
one-third of the full shot-noise value Si — lei. This re- 
sult has been credited to different theoretical approaches 
as applied to several microscopic models of disordered 
conductors. For a phase-coherent model Beenakker and 
Biittiker obtained the result using a bimodal distri- 
bution of transmission eigenvalues with the help of ran- 
dom matrix theory to calculate averages. For a semi- 
classical ID model which includes Pauli principle Na- 
gaev Pj found the same result using a Boltzmann ki- 
netic approach within an elastic and energy independent 
relaxation-time approximation. For a semiclassical se- 
quential tunneling model de Jong and Beenakker Q ob- 
tained the 1/3 value within a Boltzmann-Langevin ap- 
proach in the limit of an infinite number of equal barri- 
ers and independently from the value of their transmis- 
sion coefficient. Compatible results have been found by 
Liu et al. || from a semiclassical implementation of a 
Monte Carlo simulation which includes Pauli principle. 
For a phase-coherent model Nazarov (^] has proven the 
universality of this result in the diffusive limit for arbi- 
trary shape and resistivity distribution of the conductor 
as long as its length is greater than the carrier mean free 
path. Experimental evidence of the reduced shot-noise 
level close to the predicted 1/3 value in diffusive meso- 
scopic conductors has been provided in J^HH . 

From the above it is argued that the 1 /3 value of the 
suppression factor 7 = Si/2el is a universal phenomenon 
whose physical meaning should lay beyond classical or 
quantum mechanics and originate from some unifying 
concept. The aim of this letter is to address this issue. 
We conjecture that discreteness of charge transport is at 
the basis of such a concept, and that a transport dom- 
inated by elastic interactions is ultimately the physical 



reason for the 1/3 suppression independently from the 
quantum or classical approach used. Both the (appar- 
ently unrelated) coherent M and semiclassical |3) con- 
texts where the reduction factor 1/3 has appeared as- 
sume a degenerate Fermi gas, and the noise reduction 
comes from the regulation of electron motion by the ex- 
clusion principle. Landauer claims that the appearance 
of 1/3 in these very different cases is a numerical coin- 
cidence 0. We argue that neither phase coherence nor 
Fermi statistics are required for the occurrence of sup- 
pressed shot noise in diffusive conductors. We show a 
third case where the origin of the effect is completely 
classical, and the correlation between electrons comes 
from their Coulomb interaction, rather than the exclu- 
sion principle. The repeated occurrence of 1/3 leads us 
to believe that this must be more than a numerical coin- 
cidence. 

To support the above conjecture, we present the results 
of a Monte Carlo simulation for a nondegenerate diffusive 
conductor. By providing an exact solution of the kinetic 
equation coupled with a Poisson solver (PS) we are free 
from any approximation and thus in the position to ob- 
tain a rigorous proof of our findings. Furthermore, the 
role of long-range Coulomb interaction p~0|,[TT[ ] , although 
known since the times of vacuum diodes |12|, is here con- 
sidered for the relevant case of a medium in the presence 
of elastic and inelastic scattering for the first time. 

For the calculations we consider the following simple 
model: a lightly doped active region of a semiconduc- 
tor sample of length L sandwiched between two heavily 
doped contacts injecting carriers into the active region. 
The contacts are considered to be Ohmic (the voltage 
drop inside them being negligible) and to remain always 
at thermal equilibrium. Thus, electrons are emitted from 
the contacts according to a thermal-equilibrium Maxwell- 
Boltzmann distribution at the lattice temperature T, and 
they move inside the active region following the classical 
equations of motion by undergoing isotropic scattering 



in momentum space. To exclude additional correlations 
due to Fermi statistics, the electron gas is assumed to be 
nondegenerate. In principle, electrons are injected with 
a Poissonian statistics, i.e., the time between two consec- 
utive electron emissions is generated with a probability 
P{t) = Te~ rt , where T = ^n c VthS is the injection-rate 
density, w t h = -v/2fcsT 'firm the thermal velocity, S the 
cross-sectional area of the device, and m the electron ef- 
fective mass. However, to prove our conjecture that the 
origin of 7 = 1/3 is just elastic scattering and there- 
fore independent of the contact injection, the fluctuating 
emission rate at the contacts in the diffusive limit is taken 
to follow other kind of statistics ranging from uniform to 
Poissonian. For the simulations we have used the fol- 
lowing set of parameters: T = 300 K, m = 0.25 mo (mo 
being the free electron mass), relative dielectric constant 
e=11.7, L = 200 nm, and n c = 4 x 10 17 cm -3 (much 
higher than the sample doping). The above set of val- 
ues yields for the dimensionless parameter A = L/Luc 
(with Ld c the Debye length corresponding to n c ), which 
characterizes the importance of the electrostatic screen- 
ing |u|, the value A = 30.9, which implies significant 
space-charge effects and inhomogeneity inside the struc- 
ture. Furthermore, the average time between collisions in 
the bulk r is assumed to be independent of energy, and 
is varied from 10~ 15 s to 10~ n s, so that both regimes 
of carrier transport, ballistic {l/L ^ 1) and diffusive 
{l/L <C 1), are covered. In our calculations I, the car- 
rier mean free path, is estimated as v t hT- To analyze the 
effect of scattering inelasticity, collisions are treated as 
elastic or as inelastic. In the latter case the carrier is 
thermalized after each scattering event. 

We apply a dc voltage and calculate the time-averaged 
current / and the current autocorrelation function Cj (i) 
by means of an ensemble Monte Carlo simulator self- 
consistently coupled with a PS. We assume that the ac- 
tive region of the structure in transversal directions is 
sufficiently thick to allow a ID electrostatic treatment. 
Accordingly, the simulation is ID in real space and 3D 
in momentum space. We stress that in our approach 
the number of electrons N inside the sample fluctuates 
in time due to the random injection/extraction from the 
contacts, and we can evaluate both the time-averaged 
value (N) and its fluctuations. To analyze the impor- 
tance of the effects associated with long-range Coulomb 
interaction, we provide the results for two different simu- 
lation schemes EO ], The first one involves a dynamic PS, 
where the potential is self-consistently updated at each 
time step during the simulation by solving the Poisson 
equation under the condition that the contact potential 
remains time independent. The second scheme makes 
use of a static PS, so that only the stationary poten- 
tial profile is calculated and carriers move in the frozen 
non-fluctuating electric field profile. Both schemes give 
exactly the same average current and steady-state spatial 
distributions of all the quantities, but the noise charac- 
teristics are in general quite different. 

Figure [l] presents the results for the low- frequency sup- 
pression factor 7 = Si/2el as a function of l/L for an 
applied voltage of U = 40/csT/e. When the transport is 



ballistic or quasi-ballistic {l/L > lO^ 1 ), the static results 
show full shot noise (7^1 within numerical uncertainty) 
associated with the carrier injection (which in this range 
is modeled as Poissonian), whereas in the case of the dy- 
namic results space charge is responsible for a relevant 
noise suppression, yielding for the suppression factor the 
value of 0.045 under perfect ballistic regime {l/L S> 1) 
fllPf . The observed reduced shot-noise level is accom- 
panied by a sub-Poissonian electron number statistics 
pr| . In the transition from ballistic to diffusive regime, 
the suppression due to long-range Coulomb interaction 
remains active, being more pronounced in the case of 
inelastic scattering mechanisms. Under fully diffusive 
regime, attained at l/L < 10~ 2 , the results become in- 
dependent of the carrier injecting statistics and the fol- 
lowing asymptotic conditions are detected. The static 
results, both in the elastic and inelastic cases, keep the 
full shot-noise limit. The elastic-dynamic case attains 
the 1/3 value within numerical uncertainty in excellent 
agreement with theoretical expectations. The inelastic- 
dynamic case is further suppressed well below 1/3. 

Figure]^ complements the results of the elastic dynamic 
case in Fig. [I] by presenting calculations at increasing ap- 
plied voltages. Here we find the remarkable fact that, 
for high voltages, the 1/3 value in the diffusive limit 
is reached from the full shot-noise value in the ballistic 
limit. The reason for the different behavior in the ballis- 
tic regime is the presence or absence of a potential barrier 
near the cathode which controls the current through the 
structure. For the highest voltages, when the barrier dis- 
appears, the current saturates and the suppression fac- 
tor takes on the full shot-noisc level as for the static case 
JTo| | . It has been checked that the onset of the 1 /3 value 
takes place when, because of intensive elastic scattering, 
a significant energy redistribution among the three ve- 
locity components is achieved. Accordingly, the higher 
the applied voltage, the larger the energy gained by elec- 
trons in the free flights and, consequently, more scatter- 
ing is needed to redistribute the energy equally among 
the components. Hence, the threshold of the 1/3 regime 
is shifted towards lower values of l/L as the applied volt- 
age increases. We notice that the 1/3 limit exhibits sev- 
eral universal properties, namely, it is independent of: (i) 
the applied voltage bias (see Fig. ||), (ii) the scattering 
strength, (iii) and the carrier injecting statistics. 

To illustrate and explain the physical origin of the 1/3 
value, Fig. || reports a typical frequency spectrum of the 
suppression factor Si{f)/2el under elastic-diffusive con- 
ditions for static and dynamic PS. Monte Carlo simula- 
tion allows the calculation of the three terms in which 
Si{f) can be decomposed, namely velocity, number, and 
cross-correlation contributions Si{f) = Sy{f) + Sn{/) + 
SvN{f) jn|- In the static case the spectrum clearly 
shows that all three terms contribute to Si{f), and two 
different time scales can be identified. The longest one, 
associated with the transit time of the carriers through 
the device ry, is evidenced in the terms Sn and Svn- 
The shortest one, related to the relaxation time of elas- 
tic scattering r, is manifested in Sy- The latter is re- 
sponsible for 1/3 (within the numerical uncertainty) of 



the full shot-noise value obtained for the suppression fac- 
tor at low frequency, while the former provides the rest 
of the contribution up to the value 1. In contrast with 
this behavior, in the dynamic case the number contribu- 
tion is found to be compensated by a negative velocity- 
number cross-correlation contribution and, as a result, 
Si(f) is found to coincide with Sy(f) in all the fre- 
quency range, thus showing only the cut-off at higher 
frequencies. It is important to notice that in the fre- 
quency range 2 ^ TT ~ / ~ 5^7 both static and dynamic 
suppression factors exhibit the 1/3 value, which is related 
to velocity fluctuations. However, at low frequency only 
the dynamic case takes this value by virtue of Coulomb 
correlations, which are responsible for the mutual com- 
pensation of Sm and Svn contributions. For this com- 
pensation to take place we have checked that it is neces- 
sary to fulfil the condition L 3> Idc in order to achieve 
a significant action of long-range Coulomb interaction. 

Figure |] reports the low-frequency suppression factor 
in the diffusive regime (l/L = 10~ 3 ) as a function of the 
applied voltage, calculated by using the dynamic PS. We 
stress that these results are independent of the inject- 
ing statistics at the contacts. By comparing the elastic 
and inelastic cases we find that at the lowest bias both 
cases coincide by providing the standard thermal noise 
as predicted by Nyquist relation. On the contrary, at the 
highest voltages the elastic case achieves the 1/3 limiting 
value, in full agreement with our conjecture, while the 
inelastic case continues decreasing. 

The results for the low-frequency value of the spectral 
density in the inelastic case S} nel are closely fitted by the 
simple relation: 

S\ nel = 4k B TG j§l W 

where Go is the conductance and {N)q the average num- 
ber of electrons inside the sample, both in the limit of 
vanishing bias. This result means that, under strong 
inelastic scattering, the noise in our sample is just the 
thermal Nyquist noise (modulated by the variation in the 
number of carriers) even in the presence of a high bias 
and a net current flowing through the sample. There- 
fore, inelastic scattering strongly suppresses shot noise 
and makes the noise become macroscopic (7 « 1, see 
Fig. ^) . These results confirm previous predictions by Liu 
and Yamamoto Jl3) and Nagaev |Q , with the important 
difference that our approach is completely classical. The 
conclusion of Shimizu and Ueda Jl5| that, while dephas- 
ing is irrelevant, energy transfer from electrons to other 
systems is essential for shot-noise suppression is also sup- 
ported by present results. 

The values of the low-frequency spectral density in the 
elastic case Sf are nicely reproduced by the expression: 

which is quite similar to that obtained by Nagaev ||, 
and describes the crossover from thermal-Nyquist noise 
for eU <C k B T (Sf = 4k B TG ) to suppressed shot noise 



for eU > k B T (Sf = fel), thus providing 7=1/3 for 
the highest applied voltages (Fig. The main differ- 
ence between the effect of elastic and inelastic scattering 
on the noise in our model is the following. Strong inelas- 
tic scattering, by dissipating the energy that electrons 
gain from the field and randomizing the momentum, re- 
duces the electron average energy to that of the lattice, 
and therefore the noise corresponds to thermal-Nyquist 
noise at any bias. On the other hand, elastic scattering, 
by simply randomizing the electron momentum, allows 
for the broadening of the electron velocity distribution 
at high voltages. As a consequence, the level of noise 
increases with respect to equilibrium conditions, which 
results in a noise power with value 1/3 of the full shot- 
noise value jl6| . In addition, we have checked that if a 2D 
momentum space is considered (electron energy is shared 
between two velocity components after each scattering), 
the observed 7 = Sf/2ql in diffusive regime is about 
one-half, which indicates that the suppression factor is 
related to the momentum-space dimensionality []l7f . 

In conclusion, within an ensemble Monte Carlo scheme 
we have investigated the shot-noise suppression in non- 
degenerate diffusive conductors. Results prove that the 
1/3 value of the suppression factor is related to a trans- 
port dominated by elastic scattering (diffusive limit) un- 
der the condition eU 3> k B T. The appearance of this 
factor requires the simultaneous validity of the two con- 
ditions £ <C L and L Bc <C L, the former implying the 
achievement of fully diffusive conditions, the latter the 
decisive role of long range Coulomb correlations in sup- 
pressing the contribution associated with number fluctu- 
ations. Inelastic scattering is found to strongly suppress 
shot noise, reducing it to thermal Nyquist noise under 
heavily dissipative conditions. 
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FIG. 3. Spectrum of the shot-noise suppression factor 
for the case of static and dynamic potentials under elas- 
tic-diffusive conditions. Different contributions to the total 
value are indicated in the figure. 
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FIG. 1. Shot- noise suppression factor vs ballistic parame- 
ter IjL for an applied voltage of U = iOksT/e calculated by 
using static (open symbols, dashed line) and dynamic (full 
symbols, solid line) potentials and elastic (circles) and inelas- 
tic (squares) scattering mechanisms. 
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FIG. 2. Shot-noise suppression factor vs ballistic parame- 
ter £/L for the case of elastic scattering mechanisms at dif- 
ferent applied voltages. Calculations are performed by using 
dynamic (self-consistent) potential. 



